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Abstract 

We have optimized the parameters of extended relativistic mean-field model using a selected set of 
global observables which includes binding energies and charge radii for nuclei along several isotopic 
and isotonic chains and the iso-scalar giant monopole resonance energies for the ^"^Zr and ^^^Pb 
nuclei. The model parameters are further constrained by the available informations on the energy 
per neutron for the dilute neutron matter and bounds on the equations of state of the symmetric 
and asymmetric nuclear matter at supra- nuclear densities. Two new parameter sets BSP and 
lUFSU* are obtained, later one being the variant of recently proposed lUFSU parameter set. The 
BSP parametrization uses the contributions from the quartic order cross-coupling between oj and 
a mesons to model the high density behaviour of the equation of state instead of the oo meson self- 
coupling as in the case of lUFSU* or lUFSU. Our parameter sets yield appreciable improvements 
in the binding energy systematics and the equation of state for the dilute neutron matter. The 
importance of the quartic order uj — a cross coupling term of the extended RMF model, as often 
ignored, is realized. 

PACS numbers: 21.10.-k,21.65+f,24.30.Cz,21.60jz,26.60.+c 
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I. INTRODUCTION 



The concept of effective field tfieory has provided a modern perspective to the relativistic 
mean-field (RMF) models |l, 2|. The extended RMF models, motivated by the basic ideas 
of effective field theory, are obtained by expanding the energy density functional in powers 
of the fields for scalar-isoscalar (a), vector-isoscalar (w) and vector-isovector (p) mesons 
and their derivatives upto a given order v . The extended RMF model thus includes the 
contributions from all possible self and cross coupling interaction terms for o", u and p 
mesons in addition to the cubic and quartic self interaction terms for a meson as present 
in the conventional quantum hydrodynamic based relativistic mean field models js], 0] . The 
parameters (or the expansion coefficients) appearing in the energy density functional of the 
extended RMF model are so adjusted that the resulting set of nuclear observables agree well 
with the corresponding experimental data. The extended RMF models containing terms 
upto the quartic order [y = 4) can be satisfactorily applied to study the properties of finite 
nuclei. Inclusion of next higher order terms improves the fit to the finite nuclear properties 
only marginally. Even high density behaviour of the equation of state (EOS) of nucleonic 
matter is predominantly controlled by the quartic order 00 meson self-coupling jsj. The 
effects of higher order terms on the high density behaviour of EOS are found to be only 
modest to negligible. We would like to mention that the density dependent meson exchange 



9| and density dependent point coupling [10| versions of the RMF models are also very 
successful in describing the ground state properties of the finite nuclei. 

There have been several parameterizations {2!, [llijisl of extended RMF models. Earlier 
parameter sets Gl, G2 and TMl* were generated by considering almost all the terms upto 
the quartic order. These parameterizations yield strong linear density dependence of the 
symmetry energy coefficient and the nuclear matter incompressibility coefficient is either 
little too low or quite high. Later, FSU parameterization 16|, ll7|] demonstrated that the 
linear density dependence of the symmetry energy and the nuclear matter incompressibility 
coefficient can be reasonably obtained simply by appropriately adjusting the strengths of the 
uj — p cross-coupling and the u meson self-coupling terms of the extended RMF model. For 
the FSU parameters, the limiting mass of the neutron star is only 1.7 Mq which is somewhat 



smaller than the currently proposed limit of 2Mq 



18|. The value of the limiting mass of 



neutron stars could be significantly increased by changing uj meson self-coupling strength 
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a. 



within a reasonable range [5j. In Refs. [19|, |20|, we have shown that the strength of u 



meson self-couphng can be varied in such a way that the hmiting mass of the neutron stars 
ranges from 1.7 - 2AMq, and still bulk properties of the finite nuclei and nuclear matter at 
saturation density remain practically unaffected. 

Recently, lUFSU parameter set [2L] of the extended RMF model is obtained by readjust- 
ing the strength of the u — p cross-coupling and u meson self-coupling in such a way that 
the neutron skin of ^osp]^ nucleus is 0.16 fm and limiting neutron star mass is ~ 2.0Mq. 
The remaining parameters of the model, instead of fitting to the bulk properties of the finite 
nuclei, were tuned to yield the properties of the nuclear matter at the saturation density 
which are very much close to those for the FSU parameters. The EOSs for the symmetric 
nuclear matter (SNM) and the pure neutron matter (PNM) obtained for the lUFSU pa- 
rameters, in the density range 2.5 — 4.5po (po = 0.16 fm~^), agree reasonably with the ones 
extracted by analyzing the heavy-ion collisions data. However, energy per neutron for the 
PNM at sub nuclear densities seems somewhat larger in comparison to those obtained from 
the various microscopic approaches. The EOS for the /3-equilibrated neutron rich matter at 
high densities is softer than those deduced from the neutron star observables. The FSU or 
lUFSU like parametrizations do not include the contributions from the u — a cross-coupling 
terms. The presences of these terms might improve the high density behaviour of the EOS 
of the nucleonic matter. Nevertheless, influence of such cross-coupling on the high density 
behaviour of EOS for the nucleonic matter remain largely unexplored. 

The objective of present work is twofold. We would like to optimize the extended RMF 
model using a large set of experimental data on binding energies and charge radii instead 
of those for just a few closed shell nuclei as normally done. Our data set consists of the 
binding energies and charge radii for the nuclei along several isotopic and isotonic chains. 
The values of the constrained energies for the isoscalar giant monopole resonance (ISGMR) 
for the ^°Zr and ^°^Pb nuclei are also considered in our fit. Further, the model parameters 
are constrained by available informations on the energy per neutron for the dilute neutron 
matter and bounds on the EOSs for the SNM, PNM and the /3-equilibrated neutron rich 
matter at supra-nuclear densities. Another objective of the present work is to investigate 
whether the inclusion of u — a cross-coupling can improve the situation encountered by FSU 
or lUFSU like parameterizations. 

The paper is organized as follows. The extended RMF model is outlined briefly in Sec. Ull 
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The fit observables and the constraints employed to optimize the extended RMF model are 
discussed in Sec. IIIII In Sec. |IV]we present two newly generated parameter sets and several 
bulk properties for the symmetric nuclear matter at saturation density. In Sec. |V]we discuss 
our results for the EOSs for the nuclear and neutron star matters and compare them with 
those for few existing parameterization of the extended RMF model. The results for the bulk 
properties of finite nuclei and neutron stars obtained using our newly generated parameter 
sets are presented in Sees. IVII and [VIII Finally, summary and out look are presented in Sec. 

Ivml 



II. THE EXTENDED RMF MODEL 



The derivations of the effective Lagrangian density and corresponding ener gy d ensity 



functionals for the extended RMF model are well documented in Refs. 



12|. The 



effective Lagrangian density used in the present work can be written as j2|, |5|, 

^ = + >Co- + >Cw + + £o-ajp. (1) 

where the Lagrangian C-mm describing the interactions of the nucleons through the mesons 
is, 

I^MM = Yl ^An'd^ -(M- g,a) - {g^i^u, + ]^gp^^T.p^)]^ . (2) 

J=n,p 

Here, the sum is taken over the neutrons and protons and r are the isospin matrices. The 
Lagrangian describing self interactions for cr, a;, and p mesons can be written as, 

,2 _2\ ^3 ^ 2 3 



r ■'"/o ott 2 2\ ^3 2 3 ^^4 2 2 4 

= 2(^M^^^ -^9am^^ - ^^g^m^a , 
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24 



J^P = -Ip^.uP^" + ^rnlp^P^. 



(3) 

(4) 
(5) 



The u'^", p^^ are field tensors corresponding to the u and p mesons, and can be defined as 
^^pu _ Qfj-^i^ _ d'^u'^ and p'^'^ = d^p'^ — d^p^. The cross interactions of cr, w, and p mesons are 
described by Ca^p which can be written as. 



■'crup 



^1 

2M ^ 



4M^ 



2M 



g^rripap^p^ 



I Vlp 222 u I V2p 22 u u 

+ m^Sarripa p^p^ + j^g^m^uj^u'^p^p^. 



(6) 
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One also needs to include the contributions from the electromagnetic interaction in the case 
of finite nuclei. The Lagrangian density Cem for the electromagnetic interaction can be 
written as, 

where, A is the photon filed and F'^'^ = d^A'^ — d'^A^. The equation of motion for nucleons, 
mesons and photons can be derived from the Lagrangian density defined in Eq.([l]). The 
contributions from Eq. ([7]) are included only for the case of finite nuclei. 

It can be seen that there are five cross-coupling terms in Eq. ([6]). Two of them are the 
cubic order terms of the u — a and a — p cross-couplings and three quartic order terms 
corresponding to the u — a, a — p and u — p cross-couplings. The cross-coupling terms 
involving p-meson field enables one to vary the density dependence of the symmetry energy 
coefficient and the neutron skin thickness in heavy nuclei over a wide range without affecting 



the other properties of finite nuclei 20|, l22|, |23(]. The contribution from the u — a cross- 



couplings and self coupling of w-mesons play important role in varying the high density 
behaviour of the EOSs and also prevents instabilities in them jsl, [ill, 13 1. It may be noted 
from Eq. ([5]) that the contributions of the self-coupling of p mesons are not considered. 
Because, expectation value of the p-meson field is order of magnitude smaller than that for 



the w- meson field [2|. Thus, inclusion of the p-meson self interaction can affect the properties 
of the finite nuclei and neutron stars only very marginally 5|. 

We would like to briefiy outline the manner in which the corrections to the binding 



energies arising from the center of mass motion, pairing and quadru po 



incorporated. The spurious center of mass energy Ecm is evaluated as 24j ]. 



e correlations are 



= -17.2A-"-'MeV. (8) 

The above estimate for the Ecm is obtained by fitting the the full center of mass correc- 
tion calculated microscopically for several nuclei. We include the corrections to the binding 
energy due to the pairing correlations when the nucleon numbers are non-magic. The contri- 
butions from the pairing correlations are evaluated in the constant gap approximation with 
the gap 0, 

A = ^ MeV. (9) 
^/A ^ ' 



The pairing correlation energies for a fixed gap A is calculated using the pairing window of 
2fvjj = 82A^^^^ MeV. Soft nuclei and deformed nuclei can develop substantial contributions 



from quadrupole correlations [26|. We use a simple estimate for the quadrupole correlation 



energy or rotational correction as, 



E,ot = 2.2^/(32 - 0.05 m/m* MeV (10) 

where P2 is the nuclear quadrupole deformation and m* is the nucleon effective mass in 
bulk equilibrium matter. This estimate has been extracted by studying the microscopi cally 



computed trends of E^ot for a wide variety of Skyrme-Hartree-Fock parameterizations [27|. 
The rotation corrections are included only for the nuclei with P2 > 0.05. The Eq. fllOp is 
not relevant for the fitted nuclei as all of them are spherical. But, we shall need it for the 
binding energy systematics which is obtained by calculating the binding energies for known 
even-even nuclei and some of them are well deformed. To this end, we may point out that 
the contributions from the Coulomb exchange terms are ignored. 



III. FIT OBSERVABLES AND SOME CONSTRAINTS 



The parameters for most of the conventional and extended RMF models are obtained 



fitting the binding energies and charge rms radii only for few closed shell nuclei 



28 



i |l2, 



13 



by 



16 



291]. We fit the parameters of the extended RMF model to the experimental data for 
the binding energy, charge rms radius and energy for ISGMR. We use the binding energies 
for 62 nuclei and charge radii for 50 nuclei as listed in the Tables [T] and [Tll These nuclei lie 
along several isotopic and isotonic chains. The errors on the binding energies and charge 
rms radii used for the chi-square minimization are also given in these tables. The same set 



of nuclei are used in Ref. 30[ to obtain a set of systematically varied Skyrme forces. The 



experimental data for the ISGMR constrained energies included in our fit are 17.81 MeV for 
the ^"Zr nucleus and 14.18 MeV for the ^osp]^ nucleus 3l| with the theoretical error taken 
to be 0.2 MeV. We also constrain the model parameters, using some available informations 
on the EOS for the nuclear and neutron star matter at supra-nuclear densities together with 
energy per neutron for the dilute neutron matter. 

The EOSs for the SNM and the PNM are available in terms of the pressure versus 
density [i32j]. We considered 6 data points for each of these EOSs over the density range 
p = 2.0 — 4.5po- The EOS for the /3— equilibrated neutron rich matter are available in terms 
of pressure versus energy density jssl. We consider 32 data points for this EOS for the 
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energy density ranging from 150 — 1600 MeV/m~^. The theoretical error on the pressure 
is taken to be 25% of its required value. The realistic EOSs for the dilute neutron matter 
are available in terms of the energy per neutron versus neutron density {s^. For this case 
we take 8 data points over the range of neutron density p„ = 0.1 — 0.3po- The theoretical 
errors used for the energy per neutron is 0.05 MeV. 



IV. NEW PARAMETRIZATIONS AND NUCLEAR MATTER PROPERTIES 



We have obtained two different parameter sets of the extended RMF model. We name 
them as 'BSP' and TUFSU*', later one being the variant of recently proposed lUFSU pa- 
rameter set. The name BSP for one of our parameter set is derived from the initials of the 
authors of the present work. The parameter set BSP includes contributions from u — a 
and oj — p cross-couplings. The parameter set lUFSU* includes the contributions from u 
meson self-coupling and u — p cross-coupling. In other words, the high density behaviour of 
the EOS for the parameter set BSP is governed by the quartic order u — a cross-coupling. 
Whereas, in case of the lUFSU* it is governed by the u meson self-coupling. For the lUFSU* 
set, the strength for the u meson self-coupling is such that Co/g'i = 0.03 same as that for 
the lUFSU parameterization. This value of Co yields optimum behaviour for the EOSs of 
the SNM, PNM and neutron star matter at high densities. Similarly, for the BSP parame- 
terization, the value of 772 is adjusted to optimize the EOSs for the nuclear and neutron star 
matter at high densities. The strength ri2p of the co — p cross-coupling (Eq. (Q) is adjusted 
to yield reasonable agreement with the energy per neutron for the dilute neutron matter 
for both the BSP and lUFSU* parameter sets. The remaining parameters of the model are 
obtained by fitting the binding energies and charge rms radii of the nuclei as listed in Tables 
|T]and[ni The best fit parameters are searched using the simulated annealing method which 
has been applied to determine the parameters of the extended RMF model and Skyrme type 

In Table UTTt we list the values for the newly generated parameter 



effective forces 



19 



35 



sets. We have given the values for the Gl, G2, TMl* and NL3 parameter sets. They will 
be used to compute the EOSs for the infinite nuclear matter and the neutron star matter 
for the comparison with those obtained for the newly generated parameter sets. It must be 
pointed out that some of the coupling parameters of the parameter sets BSP, lUFSU* and 
lUFSU are quite different from the unity or equivalently they show significant deviations 
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from the naturalness. This is clearly due to the fact that not all the terms upto the quartic 
order are considered. So, the reduction in the number of the parameters is possible only at 
the cost of their naturalness. 

In Table IIVI we list the values of various quantities associated with nuclear matter calcu- 
lated at the saturation density. These quantities are evaluated as follows, 

2 (PE{p,Q) 



K = 9p 



E. 



sym 



1 £E{p,5) 
2 



d6^ 
--3p 



5=0 



dE, 



sym 



dp 



K 

^ asy 
-^sat,2 = 



Q 

2 ^ -'-'sym 



Q = 27 p- 



3 d^E{p,0) 



dp^ 



(11) 
(12) 

(13) 
(14) 
(15) 

(16) 
(17) 



where, E{p, S) is the energy per nucleon at a given density p and asymmetry S = (pn — Pp) / P- 
The incompressibility coefficient K together with -ft'sat,2 at the saturation density can yield 
the value of incompressibility coefficient for asymmetric nuclear matter 37|. One can use 
the values of -Esym, L and -K'sym at the saturation density to evaluate the density dependence 
of the s ymm etry energy coefficient which in turn can yield the EOS for asymmetric nuclear 
matter 37| . It can be seen from Table [IV] that values of -Esym and L are quite high for the 



Gl, G2, TMl* and NL3 parameterizations in comparison to those for BSP, lUFSU* and 
lUFSU parameter sets. 



V. EQUATIONS OF STATE 

We now compare our results for the EOSs for SNM, PNM and /3-equilibrated neutron 
rich matter with those obtained for a few existing parameterizations of the extended RMF 
model. As a customary we also compare our results with those for the NL3 parameter set of 
the conventional RMF model which includes non-linear terms only for the cr-mesons. It may 
be instructive to first look into the density dependence of the symmetry energy coefficient 
Egy^ and its slope L as they play predominant role in understanding the behaviour of the 

8 



EOSs for the PNM and the beta equihbrated neutron rich matter. In Figs. [T] and [2] we 
plot the density dependence of -Egym and it's slope L for various parameterizations of the 
extended RMF model as listed in Table IIIII The inset in Fig. [T] highlights the behaviour of 
-Esym at sub-nuclear densities. The Gl, G2, TMl* and NL3 parameterizations yield stiffer 
symmetry energy at supra-nuclear densities. Whereas, at densities p < 0.1 fm~^, these 
parameter sets yield softer symmetry energy. We find that -Egym = 30.0, 30.9 and 31.8 
MeV and L = 53.9, 53.9 and 49.2 MeV at po for the BSP, lUFSU* and lUFSU parameter 
sets, respectively. These values are reasonably well within -Egym = 30.5 ± 3.0 MeV and 
L = 52.5 ± 20 MeV as estimated recently by confronting the Skyrme Hartree-Fock results 
with several empirical constraints. We would like to add that the values of -Esym and L 



for the BSP, lUFSU* and lUFSU parameter sets are also in close agreement with the ones 
extracted from the experimental data on the iso- vector giant dipole resonance 0, l4o| . 

In Fig. [3] we plot the energy per neutron for the PNM at low densities. We compare 
these results with those calculated using various microscopic approaches [3^. The energy 
per neutron for the PNM for Gl, G2, TMl* and NL3 parameterizations is smaller at low 
densities in comparison to those obtained from various microscopic approaches. While, 
lUFSU yield higher values for energy per neutron for the PNM at sub-nuclear densities. 
Our parameterizations, namely, BSP and lUFSU* yield reasonable values for the energy per 
neutron at low densities for the PNM. These parametrizations include EOS corresponding to 
HE — Viowk(s'P) of Fig. [3] into the fit. The low density behaviour of the energy per neutron 
can be easily understood from the inset of Fig. [H For instance, the softer symmetry 
energy coefficient for the Gl parameter set at sub-nuclear densities is responsible for lower 
values for the energy per neutron. In Figs. H] and [5] we plot the EOSs for the SNM and 
PNM respectively in terms of pressure versus nucleon density, . The bounds on the EOSs 
as shown by shaded regions are the ones extracted by analyzing the heavy-ion collision 
data [32]. Clearly, EOSs for the BSP, lUFSU* and lUFSU parameterizations are in better 
agreement with the ones shown by shaded region. The EOS for SNM for the G2 set is also 
quite reasonable, but, it yields relatively stiffer EOS for the PNM. Thus, it appears from 
the various EOSs plotted in Figs. [3] -[5] that the overall performance of BSP, lUFSU* and 
lUFSU parameter sets are somewhat better relative to the other parameter sets considered. 

In Fig. [6|, EOSs for the beta equilibrated neutron rich matter are plotted in terms of 
pressure versus energy density. Various bounds on the EOSs as depicted by shaded regions 
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are extracted by the neutron star observables [33!]. The orange and the black boundaries of 
the shaded regions are the EOSs within la and 2a hmits, respectively. It appears at first 
glance that, except for the NL3 parameter set, all the other parameterizations yield similar 
trends for the EOS of the beta equilibrated neutron rich matter. Most of the EOSs calculated 
using extended RMF model lie within la limit for the energy densities approximately 500 — 
700 MeV/fm^. Only the EOS corresponding to the BSP parameterization stays either within 
or very close to its la limit up to very heigh energy densities. In fact, if we consider only 
the best performers, i.e., BSP, lUFSU* and lUFSU parameterizations, it can be seen that 
at very high energy densities e ^ 900 MeV/fm^, the EOSs for lUFSU* and lUFSU tend 
to go below the la limit EOS and eventually crosses 2a limit. We may recall that the 
high density behaviour for the BSP parametrization is governed by the quartic order u — a 
cross-coupling, whereas, for the lUFSU* parametrizations it is governed by the cj-meson 
self-coupling. Thus, the improvement in the high density behaviour of the EOS for the BSP 
parameter set is indicative of the importance of the contributions of the quartic order u — a 
cross-coupling. 

VI. FINITE NUCLEI 

We have computed some bulk properties of finite nuclei using the BSP, lUFSU* and 
lUFSU parameter sets of the extended RMF model. The binding energies for 513 known 
even-even nuclei are calculated to get the binding energy systematics. In Fig. [7] we display 
our results for the binding energy systematics in terms of errors 6B given as. 



where, i?th and B^xp are the theoretical and experimental values for the binding energy. We 
have indicated the present state of the art by ±1 MeV error bars. The binding energy is 
better described in case of BSP and lUFSU* in comparison to that for the lUFSU case. We 
calculate the rms error for the binding energies using the results for all the 513 even-even 
nuclei as considered in the present work. These rms errors as given in the parenthesis in the 
units of MeV are: BSP (2.7), IUFSU*(2.6) and lUFSU (5.7). It is to be noted that the rms 
error on the binding energy for our newly generated parameter sets BSP and lUFSU* are 
smaller by more than twice in comparison to that for the lUFSU. We would like to emphasize. 



SB = B, 



B, 



cxp 



(18) 



10 



this significant improvement is obviously due to the fact that we have optimized our model 
using the bulk properties of the large number of finite nuclei, unlike, the lUFSU parameters 
as discussed in Sec. [H In Ref. |27|, we presented the binding energy systematics for the 
BSR4 parameter set of the extended RMF model in which almost all the terms present in 
Eq. (E]) were included. The rms error on the binding energy for the BSR4 parameter set is 
2.6 MeV. This indicates that the terms ignored in the present work might have only small 
bearing on our results. To be more precise, the terms ignored in the present work at best 
could improve only the naturalness of the parameters. We do not show our results for the 
charge radii, since, their values for the BSP, lUFSU* and lUFSU parametrizations are very 
much similar. The rms error for the charge radii calculated using the nuclei considered in 
the fits is ~ 0.02 fm for all the three cases. 

We calculate the constraint energy Econ for the ISGMR as. 



Econ=J^ (19) 

where, mi and m^i are the energy and inverse energy weighted moments of the ISGMR 
strength function. The fully self-consistent and highly accurate values for mi and m„i are 



4l|. We follow similar strategy to 



calculated for non-relativistic mean-field models in Ref. 
evaluate mi and m_i. Once the mean-field equations are solved, the mi can be expressed 
in terms of the ground state density p as, 

mi = 2^(r2), (20) 

where, 

(r^) = / r^p{r)dr. (21) 



(22) 



The moment m_i can be evaluated via constrained RMF approach and is given as, 

where, (rj) is calculated using Eq. f l?T]) . but, the density p(r) is now obtained from the 
solutions of the constrained single-particle Hamiltonian Hx = H — Ar^.. We have compared 
our values of the i^con for ^osp]^ nucleus for the NL3 and FSU forces with those obtained 
within the RPA approach |42]. The values of i^con for both of these approaches differ at 
most by 0.04 MeV. The values of Econ for several nuclei are compared with corresponding 
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experimental data in Fig. |8l The overall trends for Econ for all the three parameter sets of 
the extended RMF model as considered are quite similar. 

The neutron skin, i?„ — Rp, the difference between the rms radii for the point neutrons 
and protons density distributions, is calculated using BSP, lUFSU* and lUFSU parameter 
sets. In Fig. [HI the values of i?„ — Rp f or the several tin isotopes are compared with the 

471]. The values of neutron-skin thickness for the 
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3^47 



corresponding experimental data 
208pb nucleus are 0.153, 0.164 and 0.161 fm for the BSP, lUFSU* and lUFSU parameter 
sets, respectively. Our values of neutron-skin for ^°*Pb nucleus are in harmony with the 
0.158!!lo'o2^fni deduced very recently from the correlations of the dipole polarizability and 
the neutron-skin 



48 



49| as well as with 0.175 ± 0.02 fm as estimated using the Skyrme 



Hartree Fock results 







VII. NEUTRON STARS 



The mass-radius relation of neutron star is important to understand the high density 
behavior of hadronic EOS. To this end, some static neutron star properties are calculated 
by solving Tolman Oppenheimer Volkov (TOV) equation. The outer crust region of the 
neutron star is described using the EOS of Riister et al. |51]. This EOS is the recent 
update of the one given by Baym, Pethick, and Sutherland [52]. Due to the fact that the 
detailed EOS of inner crust indeed is not yet certain, the polytrophic pressure-energy density 
relation is used to interpolate the EOS for the region between outer crust and the core 21 1. 
The core is assumed to be composed of either the nucleonic or the hyperonic matter in 
/3— equilibrium. The EOS of the core is obtained from the different parameter sets of the 
extended RMF model as discussed in Sec. IIVI The meson-hyperon coupling constants as 



required to compute the hyperonic EOS are taken from Refs. [20|, |53|, |5J]. In Fig. ITU] the 
mass-radius relations predicted by lUFSU, lUFSU* and BSP parameter sets are plotted 
using their nucleonic and hyperionic EOSs. For com par ison we also show the observational 
constraints extracted from the analysis of Refs. 



33 



mass observed 



18|. 



55| as well as the recent larger pulsar 



In the case of nucleonic EOS, the BSP and lUFSU* parameter sets predict slightly larger 
maximum mass compared to the ones of lUFSU. However, the BSP parameter set yields 
the radius i?i,4 for the neutron star with the canonical mass (I.4M0) which lies in between 
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those for the lUFSU* and lUFSU sets. The maximum mass predicted by each parameter 



set is in the range of the constraint region predicted by larger pulsar mass observation 



18|- 



BSP and lUFSU* predictions still also touch the upper part of 2 cr region extracted from 



the analysis of Steiner et al. 



33| . While the threshold density for the direct URCA process 



and its corresponding mass M^u as well as the transition density from core to inner crust 
and its corresponding pressure predicted by BSP and lUFSU* parameter sets are quite close 
to the ones predicted by lUFSU (see Table |V|) . 

The hyperonic star properties predicted by BSP are quite similar to the ones of lUFSU* 
except that the BSP parameter set predicts relatively smaller maximum mass radius com- 
pared to that of lUFSU*. On the other hand by comparing the particle fraction of BSP 
(right panel) and lUFSU* (left panel) of Fig. [TTl it can be seen that they have similar pat- 
tern, only the S° of BSP appears rather earlier i.e., at ~ 9 po then that of lUFSU* i.e., at 
Pb ~ 10 Pq. Thus beside difference in nonlinear term used in both parameter sets which is 
discussed in previous section, the difference in the number of TP predicted by both parameter 
sets might also influence maximum mass radius predictions. The mass-radius relationship 
as obtained using the nucleonic EOSs of the core are in better agreement with the region 
extracted from the analysis of Steiner et al. 33|- However, the maximum mass predicted 
by both parameter sets are still too small compared to the mass of pulsar J1614-2230. The 
effects of the presence of hyperons in maximum mass and Direct Urea mass predictions of 
both parameter sets can be seen in Table |Vl 



VIII. SUMMARY AND OUTLOOK 



We have optimized the extended RMF model using a large set of experimental data for 
the bulk properties of finite nuclei. The nuclear bulk properties included in the fit are the 
binding energies and charge radii of the nuclei along several isotopic and isotonic chains 
and the ISGMR energies for the ^°Zr and Pb nuclei. The density dependence of the 
symmetry energy coefficient is constrained by energy per neutron for the dilute neutron 
matter calculated in a microscopic approach. The model parameters are further constrained 
by the observational bounds on the EOS of SNM, PNM, and beta equilibrated neutron rich 
matter at supra-nuclear densities. Two different parameter sets named BSP and lUFSU* 
are obtained. Later one being the variant of recently proposed lUFSU parameter set. The 
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main difference between tlie BSP and tlie lUFSU* is tliat, in case of BSP tlie liigli density 
befiaviour of the EOS is controlled hj u — a quartic order cross-coupling and in case of the 
lUFSU* it is governed by the u meson self-coupling. We see significant improvement in the 
binding energy systematics (Fig. [7]) and the energy per neutron for the dilute neutron matter 
(Fig. [2]) over those for the lUFSU parameter set. The EOS for the beta equilibrated matter 
for the BSP parameter set (Fig. [5]) lie within or very close to its la limit bounds extracted 
using neutron star observables. These EOSs for the cases of lUFSU* and lUFSU show 
somewhat larger deviations. In particular, the pressure at higher energy densities (e > 900 
MeV/fm^) tend to become increasingly lower than the lower bound of the la limit so much 
so that it eventually crosses the 2a limit. This difference in the high density behaviour of 
the EOSs for the BSP and lUFSU* indicates importance of the contributions of the quartic 
order u — a cross-coupling, as often ignored. 

We find that the BSP and lUFSU* yield the symmetry energy coefficient Esym = 31 MeV 
and its slope L = 54 MeV and the neutron-skin for the ^"^^Pb nucleus is 0.15 — 0.16 fm. These 
values are in good agreement with the ones determined by confronting the Skyrme Hartree- 



Fock results to various empirical constraints 38|, |49|, |50|] . The neutron star properties of BSP 



and lUFSU* parameter sets using nucleonic EOS for the neutron star core are quite similar 
with the ones predicted by lUFSU. However, if hyperons are included, the BSP predict 
smaller maximum mass radius compared to the one of lUFSU*. 

The present work demonstrates that the contributions of quartic order u — a cross- 
couplings are important in order to model the high density behaviour of the EOS. The 
predictability of the extended RMF model may be still improved by including the contri- 
butions from the several gradient terms for the meson fields as presently ignored. These 
terms would contribute only for the finite nuclei allowing one to adjust simultaneously the 
properties of the finite nuclei and the infinite matter. Further, we would like to point out 
that in the present work we do not include the contributions from the 6 mesons which are 
important in order to explain the splitting of proton and neutron effective masses {56 1. The 
S meson also influences the high density behaviour of the asymmetric nuclear matter. The 
work along this lines is underway. 
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TABLE I: Experimental data for the binding energies and charge rms radii Vj-ms used in the fits 
(part I: along isotopic chain). The second line shows the globally adopted error for each observable. 
That error is multiplied for each observable by a further integer weight factor which is given in the 
column next to the data value. 
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50 


20 
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52 


20 


-436.571 








56 


28 


-483.990 


5 


3.750 


9 


58 


28 


-506.500 


5 


3.776 


5 


60 


28 


-526.842 


5 


3.818 


5 


62 


28 


-545.258 


5 


3.848 


5 


64 


28 


-561.755 


5 


3.868 


5 


68 


28 


-590.430 









18 



Table I continued. 



100 


50 


-825.800 


2 






108 


50 






4.563 


2 


112 


50 






4.596 


9 


114 


50 






4.610 


9 


116 


50 






4.626 


9 


118 


50 






4.640 


1 


120 


50 






4.652 


1 


122 


50 


-1035.530 


3 


4.663 


1 


124 


50 


-1050.000 


3 


4.674 


1 
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50 


-1063.890 


2 
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-1077.350 


2 
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-1090.400 


1 
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-1102.900 


1 
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50 


-1109.080 
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9 


5.450 


2 
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-1576.370 


9 


5.459 
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-1592.203 
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1 
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82 
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2 


5.483 


1 
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82 


-1622.340 


1 
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1 


208 


82 


-1636.446 


1 


5.504 


1 


210 


82 


-1645.567 


1 


5.523 


1 


212 


82 


-1654.525 


1 


5.542 


1 
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82 


-1663.299 


1 


5.559 


1 
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TABLE II: Experimental data for the fits, (part II: along isotonic chains). Doubly magic nuclei 
which would fit both sequences are not repeated here. For further explanations see table HI 



A 


Z 


Eb 


^rms 




±1 MeV 


±0.02 fm 


34 


14 


-283.429 


2 






36 


16 


-308.714 


2 


3.299 
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18 


-327.343 


2 
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2 
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2 
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44 


-806.849 
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-815.034 


2 
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48 


-821.064 


2 
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52 
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1 
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54 


-1141.880 


1 


4.791 


1 
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56 


-1158.300 


1 


4.834 


1 
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58 


-1172.70 


1 


4.87 


1 


142 


60 


-1185.150 
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4.915 
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62 
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2 


4.984 
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2 
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2 
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2 
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2 
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2 
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1 


210 
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1 
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-1652.510 20 
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5.555 


1 
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5.571 


1 



TABLE III: Various parameter sets for the extended RMF model. All the parameters are dimen- 
sionless. The nucleoli mass M is 939.2 MeV for the BSP and lUFSU*, 939.0 MeV for G1,G2 Q, 
NL3 [28] and lUFSU [21] and 938 MeV for TMl* Q. 





BSP 


lUFSU* 


lUFSU 


Gl 


G2 


TMl* 


NL3 


5a/4vr 


0.8764 


0.8379 


0.7935 


0.7853 


0.8352 


0.8930 


0.8131 


5a;/4vr 


1.1481 


1.0666 


1.0371 


0.9651 


1.0156 


1.1920 


1.0240 


5p/4vr 


1.0508 


0.9889 


1.0815 


0.6984 


0.7547 


0.7960 


0.7121 




1.0681 


1.1418 


1.3066 


2.2067 


3.2467 


2.5130 


1.4661 




14.9857 


1.0328 


0.1074 


-10.0900 


0.6315 


8.9700 


-5.6718 


m 


0.0872 


0.0 


0.0 


0.0706 


0.6499 


1.10 


0.0 


V2 


3.1265 


0.0 


0.0 


-0.9616 


0.1098 


0.10 


0.0 


Vp 


0.0 


0.0 


0.0 


-0.2722 


0.3901 


0.4500 


0.0 


Vip 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


V2p 
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0.0 


0.0 


0.0 


0.0 


Co 
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3.5249 


2.6416 


3.60 


0.0 


m^/M 


0.5383 


0.5430 


0.5234 


0.5396 


0.5541 


0.5450 


0.5412 




0.8333 


0.8331 


0.8333 


0.8328 


0.8328 


0.8348 


0.8333 


nip/M 


0.8200 


0.8198 


0.8216 


0.8200 


0.8200 


0.8209 


0.8126 



21 



TABLE IV: Some bulk properties of the nuclear matter at the saturation density (ps)'- binding 
energy per nucleon (B/A), incompressibility coefficient for symmetric nuclear matter {K), symme- 
try energy (£'sym)j linear density dependence of the symmetry energy (L) and various quantities 



(Ksym), (i^asy) and (i^sat2) as giveu by Eqs. dUJUSj). 



Force 


n/A 
(MeV) 


Ps 

(fm-3) 


K 

(MeV) 


77" 

(MeV) 


L 

(MeV) 


(MeV) 


J^stsy 

(MeV) 


-fVsat2 

(MeV) 


BSP 


15.9 


0.149 


230 


28.83 


50 


9 


-290 


-218 


lUFSU* 


16.1 


0.150 


236 


29.85 


50 


12 


-289 


-234 


lUFSU 


16.4 


0.155 


231 


31.30 


47 


28 


-254 


-195 


Gl 


16.1 


0.153 


215 


38.5 


123 


96 


-642 


-434 


G2 


16.1 


0.153 
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36.4 


100 


-7 


-611 


-404 


TMl* 


16.3 


0.145 


281 


37 
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-14 


-625 


-429 


NL3 


16.3 


0.148 


272 


37.4 


118 


100 


-608 


-700 



TABLE V: Tabulation of some neutron star observables in the case of without and with (+H) 
hyperons. M^^^ is the maximum mass, i?min is its radius. R\a is the radius of 1.4 Mq where Mq 
is solar mass. /3du is the threshold density for the direct URCA process, Mdu is the minimum 
neutron star's mass that may cool down by the direct URCA process, pt is the transition density 



571 ] from core to inner crust of neutron star and Pt is the pressure at pt- 





lUFSU [21] 


lUFSU* 


BSP 


IUFSU*+H 


BSP+H 




1.94 


1.96 


2.02 


1.53 


1.54 


Rrain (km) 


11.19 


11.40 


11.03 


11.65 


10.32 


RiA (km) 


12.49 


12.81 


12.64 


12.79 


12.58 


/9Du(fm"^) 


0.61 


0.61 


0.60 


0.59 


0.60 


Mbu/Mq 


1.77 


1.81 


1.82 


1.48 


1.45 


Pi(fm~3) 


0.087 


0.081 


0.087 


0.081 


0.087 


Pt(MeVfm-3) 


0.28 


0.30 


0.29 


0.30 


0.29 
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FIG. 1: (Color online) The symmetry energy coefficient Egym plotted as a function of density p/pQ 
{po = 0.16 fm^^) for several parameterizations of the extended RMF model. The inset highlights 
the behaviour of -Esvm at sub-nuclear densities. 
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FIG. 2: (Color online) Plots for the slope, L, of symmetry energy coefficient as a function of density 
for several parameterizations of the extended RMF model. 
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FIG. 3: (Color online) Plots for energy per neutron as a function of neutron density. The curves 
labeled BSP, lUFSU*, lUFSU, Gl, G2, TMl* and NL3 denote various parameter sets for the 
extended RMF model. Other curves represent the results for various microscopic approaches 34i |. 
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FIG. 4: (Color online) Pressure as a function of nucleon 



represents the EOS extracted from the analysis of Ref. 
fm^^. 



density for the SNM. The shaded area 



32( 1 ■ The density is scaled by po = 0.16 
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FIG. 5: (Color online) Same as Fig. [U but, for the PNM. 
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FIG. 6: (Color online) Pressure versus energy density for the beta equilibrated neutron rich matter. 



The shaded region represents the observational constraints taken from Ref. 



33( 1 . The orange and 



the black boundaries of the shaded regions are the EOSs within la and 2cr limits, respectively. 
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FIG. 7: (Color online) Binding energy systematics in terms of the errors 6B (Eq. [TSl) as function of 
mass number A obtained for different parameterizations of the extended RMF model. The nuclei 
that were included in the fit are marked by filled squares, well-deformed nuclei by open circles, 
and all others by triangles. Binding energy error equal to zero and ±1 MeV are indicated by 
faint horizontal lines. The corrections to the binding energies due to the pairing and quadrupole 
correlations are included for all the cases (see text for detail). 
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FIG. 8: (Color online) Fully self-consistent values for the constraint energy Econ (Eq. [T9]) of the 
iso-scalar giant monopole resonance are compared with the corresponding experimental data taken 



from Refs. 
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FIG. 9: (Color online) The neutron-skin thickness i2„ — Rp for several tin isoto pes 



lUFSU* and lUFSU parameter sets. Experimental data are taken from Ref. 



for the BSP, 



43|-l47|. 
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FIG. 10: (Color online) Mass-radius relation of BSP, lUFSU and lUFSU* parameter sets using 
nucleonic EOS and hyperionic EOS. The observational three neutron data with 1 a error bar 



that suggest small neutron star radii reported by OzellO at al. in Ref. 



55( 1 as well as two shaded 



area that suggest larger radii with 1 a and 2a error bars obtained by Steiner at al. ( 



horizontal shaded area is the mass of PSR J1614-2230 observation reported in Ref. 



33]). The 
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FIG. 11: (Color online) Particle fractions as a function of relative baryon density pb/po for lUFSU* 
and BSP parameter sets in the case hyperons are included. 
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